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1  Introduction 


Underwater  sensor  networks  consist  of  a  number  of  fixed  or  mobile  sensors,  deployed  in  the 
underwater  environment.  The  sensors  are  capable  of  sending/receiving  data  using  acoustic, 
optical,  or  radio-frequency  communication  channels  [1],  [2],  [3].  A  typical  objective  for 
such  networks  is  to  perform  data  aggregation  for  applications  as  diverse  as  environmen¬ 
tal  monitoring,  underwater  exploration,  disaster  prevention,  climate  reporting,  and  mine 
detection  [1],  [4],  [5]. 

The  acoustic  communication  channel  is  the  most  typical  physical  layer  technology  used 
in  underwater  sensors.  Unlike  the  communication  channels  used  in  the  terrestrial  sensor 
networks,  there  are  several  sources  of  uncertainty  which  influence  the  communication  be¬ 
tween  underwater  nodes  such  as  underwater  currents,  temperature  fluctuations,  multi-path 
fading,  ambient  noise,  and  sound  speed  profile  [6],  [7],  [8].  Furthermore,  these  sources  of 
uncertainty  vary  over  time  and  space  in  an  unpredictable  manner,  and  hence  one  can  use 
random  graphs  to  model  the  underwater  sensor  networks  more  effectively  [9].  This  results 
in  a  time- varying  network  structure  with  the  possibility  of  edge  addition/deletion.  Also,  the 
node  deletion  can  occur  as  a  probable  scenario  for  underwater  applications  due  to  limited 
sensor  battery  life  [10]. 

Distributed  computation  over  networks  has  attracted  considerable  attention  during  the  last 
decade  due  to  its  efficiency  in  many  cooperative  control  applications  [11].  Consensus  pro¬ 
tocols  are  among  the  fully  distributed  computation  techniques  that  are  applicable  to  data 
aggregation  tasks  in  sensor  networks  providing  robustness  to  node  failures.  In  order  to 
achieve  adequate  data  aggregation  over  a  sensor  network,  connectivity  preservation  of  the 
graph  representing  the  network  can  be  a  sufficient  and/or  necessary  condition  [12],  [13]. 
Thus,  it  is  important  to  define  a  proper  connectivity  measure  and  develop  an  efficient  algo¬ 
rithm  to  monitor,  and  if  possible  control,  the  network  connectivity  at  all  times.  Although 
there  exists  various  schemes  for  connectivity  assessment  in  the  literature  [14],  [15],  [16], 
the  corresponding  algorithms  are  not  distributed. 

An  adaptive  algorithm  for  structure  estimation  of  wireless  underwater  networks  is  devised 
in  [9].  However,  the  approach  is  constrained  merely  to  the  case  of  complete  expected 
graphs  which  imposes  some  constraints  on  the  network  communication.  This  limitation 
is  relaxed  to  certain  degree  in  the  present  publication  by  extending  the  results  to  the  case 
where  the  expected  graph  is  connected  but  is  not  complete. 

In  this  publication,  a  distributed  procedure  is  developed  to  estimate  the  connectivity  de¬ 
gree  of  an  expected  communication  graph  representing  a  sensor  network  [17].  The  sensors 
communicate  based  on  a  periodic  broadcast  schedule  in  such  a  way  that  data  interference  is 
avoided.  An  algorithm  is  employed  by  each  sensor  to  upgrade  its  estimate  of  the  expected 
communication  graph.  Moreover,  an  estimation  method  is  proposed  for  every  sensor  to 
estimate  the  probability  matrix  of  the  underlying  random  network  using  the  beliefs  of  its 
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neighbors.  Then,  weighted  vertex  connectivity  is  introduced  as  a  novel  metric  of  connec¬ 
tivity  which  extends  the  notion  of  vertex  connectivity  to  random  graphs.  The  simulation 
results  confirm  the  efficacy  of  the  proposed  strategy. 

The  remainder  of  the  publication  is  organized  as  follows.  In  Section  2,  the  problem  is 
formulated  and  some  relevant  definitions  are  given.  A  distributed  method  to  estimate  the 
expected  communication  graph  and  the  probability  matrix  of  the  network  is  introduced  in 
Section  3.  The  weighted  vertex  connectivity  degree  is  proposed  in  Section  4.  The  simula¬ 
tion  results  are  provided  in  Section  5,  and  the  concluding  remarks  are  made  in  Section  6. 

2  Problem  formulation 


The  problem  of  distributed  connectivity  assessment  for  an  underwater  acoustic  sensor  net¬ 
work  in  a  data  aggregation  application  is  investigated  in  this  publication.  Because  of  the 
random  nature  of  the  underwater  communication  channels,  some  notions  about  random 
graphs  are  borrowed  from  [12],  [13],  which  are  presented  next. 


Definition  1  Let  G  —  (V,  E)  denote  a  random  graph  composed  of  a  set  ofn  vertices  V  and 
a  set  of  directed  edges  E.  Let  also  the  matrix  P  —  \pjf  represent  the  existence  probabilities 
for  all  possible  edges,  where  plf  G  [0, 1]  is  the  probability  of  the  existence  of  the  edge 
(j,  i )  G  E.  Denote  the  adjacency  matrix  of  the  network  by  A=  [a if,  and  note  that  (j.  i)  G  E 
if  and  only  if  ap  f  0,  where  atj  is  defined  by: 

{1,  with  probability  pa, 

(1) 

0,  with  probability  1  —  Pi  j- 


Definition  2  Define  G  =  E(G)  as  the  expected  graph  of  a  random  graph  G  —  (V.E),  where 
E(-)  represents  the  expectation  operator.  The  set  of  vertices  of  the  expected  graph  G  is 
denoted  by  V  which  is  the  same  as  the  vertex  set  ofG,  and  the  set  of  its  edges  is  denoted  by 

A  A.  A. 

E.  Furthermore,  the  weighted  adjacency  matrix  ofG  is  defined  as  A  =  [djf,  where  djj  =  pij 
for  all  i,  j  G  V.  Moreover,  ( j ,  i)  G  E  if  and  only  if  ptj  f  0. 


The  sensor  network  considered  in  this  work  is  composed  of  underwater  sensors,  which  use 
acoustic  waves  for  broadcast-based  information  exchange.  In  order  to  avoid  transmission 
collisions  and  data  interference,  it  is  assumed  that  only  one  successful  broadcast  can  occur 
at  a  time.  In  other  words,  a  time  slot  is  assigned  to  each  sensor  in  a  periodic  manner 
for  broadcasting  its  data.  The  communication  graph  of  the  underwater  sensor  network 
composed  of  n  sensors  is  specified  by  a  random  digraph  G  =  {V.E),  where  its  vertex  and 
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edge  sets  are  given  by: 

V  =  {1,2 (2a) 
E  =  {(iJ)\i,jeV,aji=l}.  (2b) 

Also,  G  —  (V,E)  denotes  the  expected  communication  graph  of  the  underwater  sensor 
network,  where  V  —  V  and  its  edge  set  is  characterized  by: 

E  =  {(ij)\ijeV,Pjij=0}.  (3) 

From  the  results  given  in  [12],  a  sufficient  condition  to  reach  asymptotic  almost  sure  con¬ 
sensus  in  a  network  represented  by  a  random  graph  G  is  that  the  expected  graph  G  is 
strongly  connected.  Therefore,  data  aggregation  is  achieved  successfully  in  an  underwater 
sensor  network  as  long  as  the  expected  communication  graph  remains  strongly  connected. 
Before  introducing  appropriate  measures  for  the  connectivity  assessment  of  the  underwa¬ 
ter  sensor  network,  it  is  required  to  estimate  the  topology  of  the  expected  communication 
graph  and  its  underlying  probability  matrix  in  a  distributed  manner. 

3  Distributed  estimation  of  expected  graph 
topology  and  probability  matrix _ 

3.1  Topology  estimation 

In  the  underwater  sensor  network  under  study,  only  one  sensor  is  allowed  to  broadcast  its 
data  at  any  time  instant  in  order  to  avoid  data  interference.  A  broadcast  cycle  of  length  T 
(T  G  M>o)  is  considered,  from  which  a  time  slot  is  dedicated  to  each  sensor  to  broadcast 
its  data.  Denote  the  k- th  broadcast  cycle  by  B(k)  =  [(k  —  1  )T,  kT ),  where  k  G  N.  Let 
B'(k)  —  [(k—l)T  +  8i,  (k  —  l)T  +  8j  +  A)  be  the  k-th  broadcast  inten’al  for  sensor  i  (i  G  V), 
where  A  (A  <  E)  is  the  length  of  each  time  slot,  and  5,  G  [0,  T  —  (n  —  i  +  1)A)  indicates  the 
beginning  of  broadcasting  for  the  sensor  i  in  the  initial  broadcast  cycle  B(  \  )  =  [0.  T).  The 
acoustic  propagation  time  from  any  node  to  its  neighboring  nodes  is  taken  into  account 
by  appropriately  selecting  A.  The  broadcast  intervals  Bl(k)’ s  constitute  a  family  of  disjoint 
subsets  of  B(k)  for  every  k  G  N.  This  implies  that  \8i-8j\>A  for  all  pairs  of  distinct  nodes 

ijev. 

Each  node  broadcasts  its  estimate  of  the  global  expected  communication  graph  G  during 
its  designated  broadcast  interval.  Moreover,  the  estimate  of  G  is  updated  by  each  node 
before  its  broadcasting  starts.  The  update  procedure  for  the  /-th  node  is  performed  using 
its  previous  estimate  and  the  information  it  has  received  from  the  other  nodes  since  its 
last  broadcast.  This  time  interval  is  referred  to  as  the  k-th  receive  interval  for  node  i  and 
is  denoted  by  Rl(k).  Note  that  R’(k)  —  [{k  —  2)T  +  8/  +  A,  (k  —  l)T  +  8j )  for  all  k  >  2. 
The  initial  receive  interval  for  node  i  is  defined  as  /?'(!)  =  [0,  <5,j  for  any  i  G  V.  A  simple 
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example  is  given  in  Fig.  1  to  illustrate  the  partitioning  of  the  time  axis  for  a  network  of 
three  sensors  during  two  broadcast  cycles. 
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Figure  1:  An  example  of  periodic  broadcast  in  a  network  of  three  sensors. 


3.2  Probability  matrix  estimation 

Let  X(k)  represent  a  random  variable  with  a  Bernoulli  distribution  at  discrete  time  instant 
ke  N  such  thatX(k)  e  {0,1}.  Assume  that  po(k)  and  p\{k)  denote  the  probabilities  oftwo 
complementary  events  X  —  0  and  X  =  1  at  discrete  time  instant  k,  respectively.  Then,  the 
random  variable  X(k)  is  described  as: 


m 


1 ,  with  probability  pi(k), 
0,  with  probability  po(k) , 


(4) 


where  0  <  po(k),pi(k)  <  1  and  po(k)  +  pi(k)  =  1.  Define  X(k)  as  the  estimation  of  X(k) 
such  that: 


X(k) 


1 ,  with  probability  pi(k), 

0,  with  probability 


The  objective  is  to  design  an  estimation  procedure  such  that  the  expected  values  of  po(k) 
and  pi(k)  converge  asymptotically  to  the  real  values  of  po{k)  and  p\  (k).  respectively.  In 
other  words,  E  [pi(k)]  —>  Pi(k)  as  k—>oo  for  i  e  {0, 1}.  The  update  procedure  of  the  pro¬ 
posed  estimation  scheme  is  given  by: 


,  (k+l\  =  I  (!  -oc)pi(k)  +  a,  if  X(k)  =  i, 
PiK  }  1(1  -cc)pi(k),  if  X(k)j=i, 


(6) 


for  i  G  {0, 1},  where  a  e  (0, 1)  denotes  the  learning  rate  of  the  estimation  method.  Let 
the  probability  vector  representing  two  complementary  events  in  a  stationary  Bernoulli 
distribution  be  denoted  by  P  =  [po  p\]T .  It  is  desired  to  find  the  variance  of  the  estimate  of 
Pq  and  p\  given  in  (6),  as  time  increases. 
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Theorem  1  Consider  an  unknown  stationary  Bernoulli  distribution  given  by  (4)  and  apply 
the  procedure  (6)  to  estimate  the  probability  vector  P.  Then,  E\P(k)]  converges  to  P  as 
k  ^  co  in  an  asymptotic  manner. 

Theorem  2  Consider  an  unknown  stationary  Bernoulli  distribution  X(k)  and  apply  the 
procedure  (6)  to  identify  X(k).  Then,  the  variance  of  the  estimated  Bernoulli  distribution 
as  k—>  °o  is  given  by: 

oc 

Var[pi(° o)]  =  — -pfl  - Pi ),  (7) 

for  i  G  {0, 1},  where  a  G  (0, 1)  denotes  the  learning  rate  of  the  update  procedure. 

Consider  a  stationary  Bernoulli  distribution  with  the  estimated  probabilities  obtained  by 
using  (6),  and  let  (7)  denote  the  variance  of  the  estimated  probabilities  as  k  — >■  °o.  It  can  be 
inferred  that  as  a  — »  0,  the  variance  tends  to  zero  and  results  in  more  accurate  estimation. 
On  the  other  hand,  the  resulted  estimation  is  deteriorated  as  a  — *  1  because  of  the  increasing 
variance.  Moreover,  it  can  be  shown  that  the  learning  rate  a  G  (0, 1)  is  directly  proportional 
to  the  convergence  rate  of  E [pi(k)\  to  pt  for  i  G  {0,1}.  Therefore,  there  is  a  trade-off 
between  the  convergence  rate  and  the  estimation  accuracy  in  the  choice  of  an  appropriate 
learning  rate  a. 


Remark  1  The  results  presented  above  correspond  to  stationary  Bernoulli  distributions. 
One  can  extend  the  results  to  the  case  of  non- stationary  Bernoulli  distributions  and  derive 
an  upper  bound  on  the  estimation  error  in  terms  of  the  sampling  interval  and  the  rate  of 
change  of  the  probabilities  po{k)  and  p\  ik). 

4  Weighted  vertex  connectivity  measure 

The  objective  of  this  section  is  to  find  a  global  metric  to  evaluate  the  connectivity  degree  of 
the  expected  communication  graph  G  of  a  sensor  network.  It  is  known  that  the  strong  con¬ 
nectivity  of  G  is  a  sufficient  condition  for  asymptotic  almost  sure  convergence  to  consensus 
over  a  network  characterized  by  a  random  graph.  Consider  a  group  of  sensors  represented 
as  the  node  set  of  a  random  digraph  G  =  (V,E),  where  the  existence  probability  of  every 
directed  edge  in  E  is  assumed  to  have  a  Bernoulli  distribution.  Let  the  random  variables 
describing  the  existence  probability  of  all  edges  be  independent  and  identically  distributed. 

To  illustrate  the  importance  of  a  connectivity  measure  in  a  random  graph,  consider  the 
consensus  problem  for  two  random  networks  G\  and  G2  composed  of  three  sensors  as 
shown  in  Fig.  2. 

Let  the  purpose  of  information  exchange  between  the  nodes  in  each  graph  be  to  reach  an 
agreement  upon  a  certain  quantity  of  interest.  The  simulation  results,  as  given  in  Fig.  3, 
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Figure  2:  Two  directed  graphs  with  the  same  topology  but  different  probabilities  for  the 
existence  of  edges. 


demonstrate  faster  convergence  for  the  network  associated  with  the  random  graph  G\  which 
has  a  stronger  connectivity  due  to  higher  probability  of  the  existence  of  its  edges. 
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Figure  3:  Consensus  convergence  of  the  sensor  beliefs  for  the  graphs  G\  and  G2  in  Fig.  2. 


The  concept  of  vertex  connectivity  (VC)  has  been  used  as  a  measure  to  evaluate  the  global 
connectivity  of  digraphs  [18].  The  VC  degree  of  G  is  defined  as  the  minimum  number  of 
vertices  that  need  to  be  removed  such  that  G  is  no  longer  strongly  connected.  Let  G—(V,E) 
represent  a  digraph  with  node  set  V  and  edge  set  E.  Then  the  VC  degree  of  G,  denoted  by 
q,  is  defined  as: 

q=  nun  q(i,j),  (8) 

t,j£V,  1 


where, 


N(i,j),  if  (iJ)£E, 
\V\-h  if  (iJ)eE, 


(9) 


and  N(i,  j )  denotes  the  maximum  number  of  vertex-disjoint  directed  paths  connecting  i  to 
j  in  G. 
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Several  polynomial-time  algorithms  are  provided  in  the  literature  to  find  the  VC  degree  of  a 
graph  [19],  [20],  [21].  The  general  idea  behind  these  algorithms  is  that  the  minimum  num¬ 
ber  of  vertices  whose  removal  disconnects  any  pair  of  non-adjacent  vertices  is  equal  to  the 
maximum  number  of  mutually  vertex-disjoint  directed  paths  between  them  (see  Menger’s 
Theorem  [18]). 

However,  this  measure  does  not  account  for  the  probability  matrix  P  and  this  calls  for  a 
more  accurate  measure  of  connectivity  capturing  the  probabilistic  nature  of  the  environ¬ 
ment.  The  weighted  vertex  connectivity  (WVC)  measure  extends  the  notion  of  the  VC 
degree  to  the  more  general  case  of  weighted  graphs.  This  nonnegative  measure  is  positive 
for  a  strongly  connected  graph,  and  a  larger  value  of  this  measure  represents  “stronger” 
connectivity.  To  clarify  this  new  concept,  the  multiplicative  path  weight  is  defined,  which 
is  based  on  the  mutual  independence  of  the  random  variables  used  for  describing  the  prob¬ 
abilistic  nature  of  the  network’s  edges. 


Definition  3  Let  the  ordered  set  of  distinct  nodes  { z'o .  i\ .  denote  a  directed  path  from 

z'o  to  im  in  the  expected  communication  graph  G.  The  multiplicative  path  weight  of  this 
directed  path,  denoted  by  PW  (z'o,  im),  is  defined  as  follows: 

m 

PW  (z'o,  im )  =  n  Pilch-  v  ( 1  °) 

k=  1 

where  ptj  is  the  (z,  j)  element  of  the  probability  matrix  P.  Furthermore,  the  length  of  the 
above  path  is  given  by: 

PL(i0,im)=m.  (11) 


Since  every  element  of  the  matrix  P  represents  the  probability  of  the  existence  of  the  corre¬ 
sponding  edge  in  the  expected  communication  graph  G,  the  multiplicative  path  weight  can 
be  interpreted  as  the  probability  of  the  existence  of  a  given  path,  as  long  as  the  probabilities 
of  the  existence  of  different  edges  of  the  path  are  mutually  independent.  In  order  to  extend 
the  notion  of  vertex  connectivity  degree  to  random  graphs,  it  is  first  required  to  define  the 
local  WVC  measure  for  any  pair  of  distinct  nodes  z,  j  €  V  in  the  expected  communication 
graph  G  =  (V,E).  This  local  measure,  denoted  by  q(i,j),  is  defined  as  the  maximum  value 
of  the  summation  of  the  multiplicative  weights  of  all  vertex-disjoint  directed  paths  from  z 
to  j  which  belong  to  G.  In  other  words,  q(i.  j)  represents  the  maximum  of  the  summation 
of  the  existence  probability  of  the  vertex-disjoint  paths  from  i  to  j.  Let  T(z,  j)  contain  all 
possible  directed  paths  from  i  to  j  in  G.  Then: 


q(ij) 


N(iJ ), 

max  ((|V 


if  (z,  j)  f  E, 
l)pij,N(i,j)),  if  (i,f)  G  E, 


(12) 
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where, 


N(iJ)=  max  £  PWk(iJ ),  (13) 

and  S(i,  j)  represents  a  set  of  directed  paths  from  i  to  j  which  are  mutually  vertex-disjoint. 
Then,  the  WVC  metric,  denoted  by  q ,  is  introduced  as  a  global  connectivity  measure  of  the 
expected  communication  graph  G  as  follows: 

q=  min  q{i,j).  (14) 

i,j€V,ij=j 

The  WVC  measure  q  can  be  considered  as  an  extension  of  the  VC  degree  q  for  weighted 
graphs,  where  an  edge  weight  in  the  present  problem  denotes  the  existence  probability  of 
that  edge,  and  varies  between  zero  and  one.  Furthermore,  the  WVC  measure  q  is  a  non¬ 
negative  real  value  ( q  G  R>o),  and  is  more  sensitive  to  changes  in  the  network  compared 
to  the  VC  degree  q  which  is  a  nonnegative  integer  (<7  €  Z>q).  Note  that  For  two  differ¬ 
ent  expected  communication  graphs  G\  and  G2  with  <71  >  qi ,  it  is  possible  that  <71  <  <72, 
depending  on  the  probability  matrices  P\  and  IP.  Note  also  that  since  the  random  nature 
of  the  environment  is  captured  by  the  WVC  measure,  it  is  more  suitable  for  the  networks 
which  are  represented  by  random  graphs. 

4.1  Illustrative  examples  to  compare  the  connectivity 
measures 

Example  1:  As  the  first  example,  consider  the  expected  communication  graph  G  shown 
in  Fig.  4,  where  the  elements  of  the  probability  matrix  are  shown  as  weights  on  all  edges. 
Using  the  definition  of  VC  degree  and  on  noting  that  G  is  strongly  connected,  one  obtains 
<7=1.  Based  on  the  definition  of  WVC  measure,  q(a,c)  —  q(c,a)  =  0.72,  which  results  in 
<7  =  0.72. 


Figure  4:  Expected  communication  graph  of  Example  1 

Example  2:  In  the  second  example,  the  expected  communication  graph  G  shown  in  Fig.  5 
is  considered.  This  graph  is  similar  to  the  one  used  in  Example  1  with  the  difference  that 
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the  values  of  edge  weights  pag-,Pga-,PbciPcb  have  been  changed  from  0.4  to  0.3  in  this 
example.  Like  the  previous  example,  the  edge  weights  shown  in  this  figure  represent  the 
existence  probability  of  the  corresponding  edges.  Since  the  graph  is  strongly  connected 
and  there  exists  at  least  one  directed  path  between  any  pair  of  its  distinct  nodes,  G  is  1- 
vertex  connected  (q  =  1)  similar  to  Example  1.  Unlike  Example  1,  the  minimum  local 
WVC  measure  is  given  by  q{a,c)  =  q(c,a)  =  0.648,  which  leads  to  q  =  0.648. 


Figure  5:  Expected  communication  graph  of  Example  2 


Remark  2  The  results  developed  in  this  work  can  be  used  to  improve  the  performance  of 
an  underwater  acoustic  sensor  network  by  identifying  the  relative  contribution  of  each  link 
in  the  connectivity  of  the  ovemll  network.  For  example,  consider  an  underwater  sensor 
network  represented  by  the  random  graph  given  in  Fig.  6.  While  there  are  two  edges  of 


Figure  6:  Expected  communication  graph  of  Remark  2 

the  same  probability  0. 1  in  this  graph,  it  results  from  the  proposed  connectivity  measure 
that  the  edge  from  node  3  to  node  1  plays  a  more  important  role  than  the  edge  from  node 
1  to  node  3  in  the  connectivity  of  this  graph  because  there  exists  another  path  from  node 
1  to  node  3  with  maximum  probability  of  existence  while  the  edge  from  node  3  to  node  1 
represents  the  only  path  from  node  3  to  node  1. 
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5  Simulation  results 


Consider  a  network  of  six  underwater  acoustic  sensors  with  the  periodic  broadcast  schedule 
described  in  Section  3.  Assume  that  the  existence  probability  of  the  edges  in  the  commu¬ 
nication  digraph  G  is  described  by  a  time-varying  matrix  P(k)  given  below: 

'  0  0.5  0  0.8  0.7  0 

0.6  0  0.6  0  0.5  0.9 

p/i\ _  0.7  0.6  0  0.9  0  0 

0  0.7  0.4  0  0.9  0.7 

0  0.7  0  0.9  0  0.7  +  0.3  sin  cok 

_  0.8  0  0.3  0.9  0  0 

where  (0  —  0.005,  and  k  e  N  denotes  the  k-th  broadcast  cycle  (note  that  the  above  matrix  is 
only  used  in  simulations,  and  is  considered  as  an  unknown  matrix  in  all  algorithms).  The 
expected  communication  graph  G,  which  is  constructed  based  on  the  random  graph  G  with 
the  edge  weight  matrix  P(k),  has  the  following  adjacency  matrix: 

0  10  110 
10  10  11 
110  10  0 
0  110  11 
0  10  10  1 
10  110  0 

The  expected  communication  graph  G  induced  by  A  is  depicted  in  Fig.  7. 


Figure  7:  Expected  communication  graph  G. 

Let  the  length  of  each  broadcast  cycle  be  T  —  6  sec,  and  the  length  of  each  broadcast 
interval  be  A  =  1  sec.  Let  also  the  broadcast  order  of  the  i-th  sensor  be  specified  by  <5,  = 
i  —  1  for  i  G  (1,2,  ...,6}.  Using  a  —  0.01  as  the  learning  rate,  the  estimated  network  size 
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perceived  by  each  node  versus  the  number  of  broadcast  intervals  is  shown  in  Fig.  8.  On 
the  other  hand,  Figs.  9  and  10  present  estimates  of  the  edge-set  size  and  the  VC  degree, 
respectively,  by  each  node.  Fig.  8  demonstrates  that  all  nodes  identify  the  network  size  after 
3  broadcast  cycles,  while  Figs.  9  and  10  show  that  all  nodes  precisely  estimate  the  edge-set 
size  and  VC  degree  after  10  broadcast  cycles  (note  that  each  broadcast  cycle  consists  of  six 
broadcast  intervals). 


Figure  8:  Network  size  estimation  perceived  by  every  sensor  versus  the  number  of  broad¬ 
cast  intervals. 

In  order  to  verify  the  effectiveness  of  the  probability  estimation  method  and  the  proposed 
WVC  measure,  the  corresponding  algorithms  are  simulated  for  a  sequence  of  1000  broad¬ 
cast  cycles.  The  results  of  the  probability  estimation  procedure  to  identify  the  constant 
element  p  15  and  the  time-varying  element  p^  are  demonstrated  in  Figs.  11  and  12,  respec¬ 
tively.  Fig.  11  shows  the  efficacy  of  the  estimation  procedure  for  a  stationary  Bernoulli 
distribution,  while  Fig.  12  shows  the  performance  of  the  proposed  method  in  estimating 
a  non- stationary  Bernoulli  distribution.  In  Fig.  13,  the  connectivity  measures  q  from  the 
viewpoint  of  sensor  1  is  demonstrated. 

6  Summary 


Connectivity  assessment  of  an  underwater  sensor  network  using  distributed  algorithms  is 
investigated  in  this  publication.  The  sensors  are  assumed  to  broadcast  their  information  in  a 
periodic  manner  using  acoustic  modems,  and  the  underlying  communication  graph  is  mod¬ 
eled  by  a  random  graph  due  to  unpredictable  conditions  of  the  underwater  environment. 
An  update  procedure  is  used  by  sensors  to  identify  the  structure  of  the  expected  commu¬ 
nication  graph  at  first.  Also,  a  novel  estimation  procedure  is  introduced  to  estimate  the 
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Number  of  Broadcast  Intervals 


Figure  9:  Edge-set  size  estimation  perceived  by  every  sensor  versus  the  number  of  broad¬ 
cast  intervals. 

probability  matrix  in  a  distributed  manner  by  each  sensor.  Then,  the  weighted  vertex  con¬ 
nectivity  is  proposed  as  a  novel  measure  to  evaluate  the  connectivity  of  the  network  which, 
in  fact,  extends  the  vertex  connectivity  metric  to  the  case  of  random  graphs,  capturing  the 
probabilistic  nature  of  the  underwater  sensor  networks.  The  main  challenge  for  the  future 
work  is  to  come  up  with  computationally  efficient  algorithms  to  assess  the  connectivity 
of  a  random  network.  Simulation  results  demonstrate  the  effectiveness  of  the  algorithms 
developed  in  this  publication. 
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Figure  10:  VC  degree  perceived  by  every  sensor  versus  the  number  of  broadcast  intervals. 
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Figure  11:  Estimated  and  actual  probability  perceived  by  every  sensor  versus  the  num¬ 
ber  of  broadcast  cycles. 
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Figure  12:  Estimated  and  actual  probability  perceived  by  every  sensor  versus  the  num¬ 
ber  of  broadcast  cycles. 


Figure  13:  WVC  metric  q  from  the  viewpoint  of  sensor  1. 
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